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Measurements of Electroweak Penguin Processes b —^ d,sj and b — > s£+£ at Belle 
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We report the measurements of electroweak penguin processes b ^>- d,sj and b -^ s£^£~ at 
the Belle experiment. The analyses are based on 152 million BB events collected at the T(46') 
resonance with the Belle detector at the KEKB e^e~ asymmetric-energy collider. The b ^ dj 
transition is studied through exclusive decays: B" -^ p^y, S" — * lo^ and B'^ -^ P^7- We 
apply fully inclusive and semi-inclusive reconstruction to 6 -^ S7 and b -^ s£^£~ processes, 
respectively. The results on the exclusive b -^ d'y and semi-inclusive b -^ st^l~ measurements 
are preliminary. 



1 Introduction 



In the Standard Model (SM), flavor-changing neutral current (FCNC) process is forbidden at 
tree level. Such a process can occur at higher order via electroweak penguin and box diagrams. 
The FCNC process is therefore sensitive to the effects of non-SM physics that may enter these 
loops. 

In this report, we show the measurements of FCNC processes b ^ d,s^ and b -^ si~^l~. 
The b ^>- dj process is studied through exclusive decays: B^ -^ p^j, B^ — > uj^ and B^ — > p'^'f. 
We apply fully inclusive and semi-inclusive reconstruction to 6 ^ 57 and b — > s£"^^~ processes, 
respectively. These measurements are based on a sample of 152 x 10^ T(4S') -^ BB decays 
collected with the Belle detector^ at the KEKB^ e'*"e~ asymmetric-energy collider. 



2 Exclusive b ^ dj 



The b —I- dj process is suppressed with respect to 6 — > 57 by the Cabibbo-Kobayashi-Maskawa 
(CKM) factor |Vid/Visp ~ 0.04 with a large uncertainty due to lack of precise knowledge on Vtd- 



In this analysis, we use the exclusive reconstruction of the decays, B~^ 



P"7: 



5° 



P°7, 



and B^ 



UJJ. 



We also reconstruct B~^ -^ K*~^^ and B^ — > K*^^ for control samples. The 



following decay chains are used to reconstruct the intermediate states: /9^ 



TT+vr*^, p^ —> vr+TT 



UJ 



TT TT TT 



K*+ -^ K+7r°, K*^ -^ K+TT', and vr" 



77. In each event, we select the photon 
with the largest energy in the range 1.8 GeV < E^ < 3.4 GeV, in the e^e~ center-of-mass 
frame (CM). To suppress backgrounds from vr*^ -^ 77 and rj ^ '^^ decays, we veto the event 
if the reconstructed mass of the primary photon and any other photon of 30 (200) MeV or 
more is within ±18 (32) MeV/c^ of the ir^ {rf) mass. B candidates are formed by combining 
a /9 or a; candidate and the primary photon using two variables: the beam-energy constrained 

|p^/cp and the energy difference /S.E 



mass Mbc = ^{E^^^jW 



E 



B 



^bcam' where 



p*^ and E'^ are the measured CM momentum and energy, respectively, of the B candidate, 
and -Ebeam is ^^^ C^ beam energy. The photon energy is replaced by -Ebcam ~ ^*p/uj ^^ ^^^ 
momentum p*^ is calculated. We define the signal region as —0.1 GeV < AE' < 0.08 GeV and 
5.273 GeV/c^ < Mbc < 5.285 GeV/c^. 

There are two major background sources from B decays: B -^ K*j, and B -^ p/co7r^. 
To suppress B -^ K*j, we calculate Mk-k, where a kaon mass is assigned to one of the pion 
candidate, and reject the candidate if Mk-k < 0.96 (0.92) GeV/c^ for the p^^ {p^l) mode. To 
reject B -^ p/ujn^, we apply a helicity angle cut |cos0hcil > 0-8 (0.6) for p^'^ and 0^7 (/9^7) 
modes. Here, ^hei is the angle between tt"*" and B momentum vectors in the p rest frame, or 
between the normal to the u decay plane and the B momentum vector in the to rest frame. 

The background from continuum e+e~ —^qq{q = u, d, s, c) events is rejected using the event 
topology information. We use a Fisher discriminant? constructed from 16 modified Fox- Wolfram 
moments ^ and the scalar sum of the transverse momentum ^ . We also use the decay vertex of 
the candidate B meson as well as the origin of the remaining tracks in the event. The difference 
between these two vertices along the z-axis discriminates continuum events that has a common 
decay vertex and signal events whose decay vertices are displaced in the laboratory frame. 

To obtain the signal yield, we perform an unbinned maximum likelihood fit to Mbc and AE 
to the events in \A.E\ < 0.3 GeV and Mbc > 5.2 GeV/c^. We perform a simultaneous fit to 
three signal modes (B -^ (p, uj)j) plus two B -^ K*^ modes, assuming the isospin relations 
B{B -^ (p,w)7) = B{B+ -^ p+-f) = 2^B{B^ -^ pS) = 2^B{B^ -^ ujj) and B{B+ -^ 

K*+7) = l£±J3{B'^ -^ K*^l), where we use ^ = 1.083 ± 0.017^. Here, floated parameters are 
the branching fractions for B -^ (p, lo)^ and B -^ K*^, five background normalizations and five 
background AE slopes. 

Preliminary results of the simultaneous fit are shown in Fig. 1 and given in Table 1. We find 
the simultaneous fit gives a significance of 3.5, where significance is defined as y^ —2ln{Co / Cmax) , 
-^max is the maximum likelihood in the Mbc fit, and Cq is the likelihood of the best fit here the 
signal yield is constrained to be zero. 



'Throughout the article charge-conjugate states are imphcitly included. 



Mode 



efficiency 

(isyst) 



signal yield 

(ibstatibsyst) 



significance 



branching fraction 

(±stat±syst) 



(1.8+^;^ ±0.1) X 10-^ 
(0.5t{]i±0.2) X 10^^ 
(1.3t{]-]±0.2) X 10^6 
(1.8;gj±0.1) X IQ-^ 



B+ 


-> p+-f 


(5.6 ±0.4)% 


15.5+^:^ ±1.5 


2.5 


B°- 


-A 


(5.0 ±0.3)% 


q f? +3.6 +0.7 
•-*-"-2.8 -0.9 


1.2 


B^- 


-^ ix'7 


(4.7 ±0.5)% 


8.9ltg±1.2 


2.3 


B- 


> {p,uj)-f 


— 





3.5 



Table 1: Results of the efficiency, signal yield, significance, and branching fraction from simultaneous and indi- 
vidual fits. These numbers are preliminary. 




Figure 1: Projections of the simultaneous fit results to Mhc and A_E of the individual modes to the signal region, 

for B — > LJ7 (a and b), B^ — > p^7 (c and d), and B -^ p 7 (e and f). Lines represent the total fit result (black), 

signal (red), K*^ (blue) and continuum (green) components. B — » (p/uj)-K^^ components are invisibly small. 



3 Inclusive b 



sj 



We recently measured the b ^ s'j branching fraction using a fully-inclusive approach. The 
detailed descriptions are written in Ref. ''' In this analysis, to extract the signal 6 — > ,37 spectrum, 
we collect all high-energy photons, vetoing those originating from vr'^ and t] decays to two photons. 
The contribution from continuum events is subtracted using the off-resonance sample. The 
remaining backgrounds from BB events are subtracted using Monte-Carlo (MC) distributions 
scaled by data control samples. After subtracting the backgrounds, we correct the photon 
spectrum for the signal selection efficiency function obtained from signal MC, applying the 
correction determined by data control samples. 

The efficiency-corrected spectrum is shown as a function of CM photon energy in Figure 2. 
The two error bars for each point show the statistical and the total error, including the sys- 
tematic error which is correlated among the points. As expected, the spectrum above the 
3 GeV endpoint for decays of B mesons from the T(4S') is consistent with zero. Integrating 
this spectrum from 1.8 to 2.8 GeV, we obtain a partial branching fraction of B{b — > S7) = 
(3.59 ± 0.32 1q'3^ -oio?) ^ 10~^, where the errors are statistical, systematic and theoretical, re- 
spectively. This result is in good agreement with the latest theoretical calculations^'^. We have 
also measured the moments of the distribution and obtain {E^) = 2.289 =b 0.026 =b 0.034 GeV 

and (^E^'j - {E^f = 0.0311 ± 0.0073 ± 0.0063 GeV^ for E* > 1.8 GeV, where the errors are 
statistical and systematic. 



25000 




-5000 



E\ [GeV] 



Figure 2: Efficiency-corrected photon energy spectrum. The two error bars show the statistical and total errors. 



4 Semi-inclusive b 



d+e- 



We study the inclusive B -^ Xgi^i" process, where i is an electron or a muon and Xg is a 
hadronic system containing an s-quark. The s-quark hadronic system is reconstructed with a 
semi-inclusive reconstruction approach: reconstruct with one K^ or Kg and up to four pions 
(at most one pion can be neutral). Compared to a fully inclusive approach, this method has the 
advantage of having the strong kinematical discrimination by using the beam-energy constrained 
mass Mbc and the energy difference AE. 

In addition to the discrimination, further background suppression to reduce the large com- 
binatorial backgrounds is necessary. The main contribution to the combinatorial background 
is from semileptonic decays in BB events. In these events, B — > Xsi~^i~ candidates are re- 
constructed with the decay products from both BB mesons . Such a background event has a 
significant amount of missing energy due to the neutrinos from the semileptonic decays. Another 
contribution to the combinatorial background is from continuum events, which are effectively 
suppressed with event-shape variables. 

There are two background sources that can peak in Mbc and AE. The first is from 
B -^ J/ipX and B -^ il}{2S)X decays with J/ip{tp{2S)) -^ £^£~. This peaking background 
is efficiently removed with cuts on the dilepton mass m{l^i^). The resulting veto sample pro- 
vides a large control sample of decays with a signature identical to that of the signal. The 
second is from B -^ K^{Kg)mT (n > 1) decays with misidentification of two charged pions as 
leptons. We estimate these peaking background contaminations, then subtract them from the 
signal yield. 

We perform an extended, unbinned maximum likelihood fit to the M^c distribution in the 
region M^c > 5.2 GeV/c^ to extract the signal yield as well as the shape and yield of the 
combinatorial background. We fit the Mi,c distributions for the selected B -^ Xse^e~ and 
B —> XsH^fi~ candidates separately, and B — > Xsi~^i~ {i = e,fj,). The obtained results are 

XsC fi^ candidates. 



10^ BB decays, we 



shown in Figures 3(a)-(c). Figure 3(d) shows the distribution for B - 
There is no peaking component as expected. Using a sample of 152 
observe a signal of 72.3 ± 13.8{stat) ± A.6{syst) events. 

Figures 4(a)-(c) show the distributions of hadronic mass M^s, lepton pair mass Mu, and q^ = 
M^ for electron and muon channels combined, obtained by performing the nominal likelihood 
fit in separate M^s, Mu, and q^ regions. 



(a) Xse+e- 
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Figure 3: Distributions of Mf,c for selected (a) B -^ Xse^e~ , (b) B -^ Xs^i^ fj.~ , (c) B —^ Xs£'^i~ , {(. — e,n), and 

(d) B — > XgC'^fj,^ candidates. The red lines represent the result of the fits, and the green and blue lines represent 

the peaking and combinatorial background components under the signal peaks, respectively. 
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Figure 4: Distributions of number of signal events as a function of (a) hadronic mass Ai^s (b) lepton pair mass Mu, 
and (c) q^ = M^ for electron and muon channels combined for data (points) and Monte Carlo signal (histogram). 

The vertical error bars represent statistical errors only. 



The preliminary results of the measured branching fractions for M[i'^l ) > 0.2 GeV/c^ are 

B{B -^ XsS+e-) = (4.45 ± l.S2tojt) x IQ-^, 
B{B -^ X,/i+/i^) = (4.30 ± imtllf) X 10"^ and 
B{B -^ Xsl+e-) = (4.39 ± 0.84+0;^^) x 10-^ 

where the first error is statistical and the second error is systematic. The combined B -^ 
Xsi~^l~ branching fraction is the weighted average of the branching fractions for the electron and 
muon channels, where we assume the individual branching fractions to be equal for M{i'^i~) > 
0.2 GeV/c2. 

5 Summary 

We have studied the electroweak penguin processes b ^ d,sj and b — > si^i^ based on a sample 
of 152 X 10^ T(4S') —f BB decays at the Belle experiment. The b ^ d'y transition is studied 
through exclusive decays: B^ -^ p^^, B^ -^ loj and B^ -^ P^l- We apply fully inclusive and 
semi-inclusive reconstruction to 5 ^ 57 and b -^ si^i^ processes, respectively. 

We present the first evidence for the b ^ d'j process using a simultaneous fit to the B ^ pj 
and B^ -^ oj^ modes. The preliminary result is consistent with the SM predictions lO-ii. We 
measure the branching fraction and photon energy spectrum of inclusive b ^ S'j process in the 
energy range 1.8 GeV < E* < 2.8 GeV. For the first time 95% or more of the spectrum is 
measured, allowing the theoretical uncertainties to be reduced to a very low level. We measure 
the branching fraction for the B -^ Xsi'^i~ , where i = e or p and Xg hadronic system is semi- 
inclusively reconstructed. The preliminary measurement is consistent with the recent prediction 
by Ah et al.^^. 
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